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We investigate theoretically an open dynamics for two modes of electromagnetic field inside a 
microwave cavity. The dynamics is Markovian and determined by two types of reservoirs: the 
"natural" reservoirs due to dissipation and temperature of the cavity, and an engineered one, provided 
by a stream of atoms passing trough the cavity, as devised by Pielawa et al. [Phys. Rev. Lett. 98, 
240401 (2007)]. We found that, depending on the reservoir parameters, the system can have distinct 
"phases" for the asymptotic entanglement dynamics: it can disentangle at finite time or it can have 
persistent entanglement for large times, with the transition between them characterized by the 
possibility of asymptotical disentanglement. Incidentally, we also discuss the effects of dissipation 
on the scheme proposed in the above reference for generation of entangled states. 
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I. INTRODUCTION 

In recent years, the problem of entanglement dynam- 
ics has gained the attention of the quantum-information 
community jlj. Despite the fact that some recent re- 
sults contradict the intuition that "the more entangled 
the better" Q, to protect entanglement is yet an impor- 
tant concern when one wants to make quantum comput- 
ers. The first works on the "fate" of entanglement on 
an open system focused on two or more qubit systems. 
Curious names such as entanglement sudden death have 
appeared, first by the recognition of the possibility of en- 
tanglement between two qubits to vanish in finite time, 
although coherences only decay exponentially. Similar 
phenomena also happen for continuous variables (CV) 
systems A general geometrical picture has already 
been offered showing that the long term behavior of en- 
tanglement depends essentially on the set of asymptotic 
states 0| . 

In a series of papers, Paz and Roncaglia studied some 
phase diagrams for the asymptotic behavior of entan- 
glement on two-mode Gaussian states (GS) exposed to 
a common thermal environment Q. In particular, they 
built them in terms of squeezing and temperature for two 
harmonic oscillators. For the resonant case, they found 
three very distinct possible fates: entanglement can suf- 
fer from sudden death, can enter in a perpetual cycle 
of death and birth, or can be persistent. From a differ- 
ent perspective, Pielawa et al. made a proposal of how 
to use an engineered reservoir to create two-mode Gaus- 
sian entangled states in cavity quantum electrodynamics 
(CQED) @. 

In this paper we make a threefold study: we gener- 



alize the studies from Paz and Roncaglia, by construct- 
ing phase diagrams in terms of different variables, for 
systems subjected to not only natural reservoirs, but 
also this engineered one; we show how those phase di- 
agrams can be experimentally obtained: and incidentally 
we study the robustness, against thermal noise, of the en- 
gineered reservoir strategy for generating two-mode en- 
tangled states. 

The next section is devoted to reviewing some basics 
about two-mode GS and the proposal for generating en- 
tangled ones in CQED @. Section ITU is the central part 
of our study, where the entanglement dynamics of two 
modes subjected to thermal noise and the engineered 
reservoir is discussed. The experimental proposal is fo- 
cused in Sec. II VI followed by the study of the robustness 
of the method suggested in Ref. (H . Discussion and con- 
cluding remarks close the paper. 



II. PRELIMINARIES 

Here we will review some necessary definitions and 
tools in the theory of entanglement of GS and also the 
scheme to implement the engineered reservoir. 

A. Entanglement in Gaussian States 

Many of the protocols and techniques from the en- 
tanglement theory of finite dimensional systems can be 
adapted to CV systems, usually restricted to the set of 
GS. For instance, the protocols of teleportation 0, HI 
and quantum key distribution 0, [l^ have their analogs 
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for CV systems (and may be more robust experimentally 
flll|); the Peres-Horodecki fl2l , f]~3l | separability criteria 
can be applied, and is also necessary and sufficient for 
two- mode GS fl4) ; and entanglement quantifiers such as 
negativity fl6| . logarithmic negativity fl7| , and entan- 
glement of formation [l8j can be computed within the 
formalism of symplectic geometry |l9l42l| . 

GS are defined as those states whose Wigner charac- 
teristic function is Gaussian, so they are completely de- 
scribed by its first and second statistical momenta. First 
momenta, mean values, can be locally changed by re- 
defining the mode operators, so all entanglement infor- 
mation is given by second momenta. Choosing only two 
modes, with destruction operators a,j = (ii + ipi)/V2, 
i = 1,2, and corresponding quadrature amplitudes Xi,pi, 
the Wigner characteristic function of a state p is given by 
x(z\,Z2) = Tr[pcxp (zi&i — z^a\ + Z2&2 — z^al)], where 
Zi are complex numbers. The state second momenta, on 
the other hand, are well grouped under its covariance 
matrix (CM): 

f ni + 5 mi m s m c \ 

1 m s m c ri2 + 5 m 2 

\ m* m s m* 2 n-x + \ ) 

where rij = (ajdj), rrij = — (a 2 ), m s = — (a\a\), m c = 
(0.102), and (£) denotes the quantum expectation value 
Tr£/5 of an observable £. The state is Gaussian (with null 
first momenta) if and only if x( z ii z 2) = exp (—^z^VpZ). 
where z* = (z\, z*, Z2, z%)- 

We can write the CV as a block matrix, as follows: 

where Vi is a 2 x 2 matrix related to the mode i, and 
C is a 2 x 2 matrix that gives the correlations (both 
quantum and classical) between the modes. One should 
note that to represent a quantum state, a CM should also 
obey the generalized Robertson-Schrodinger uncertainty 
relations [21|. A CM not obeying such relations is called 
nonphysical. 

Simon [l4| has shown that, as for two qubits, the Peres- 
Horodecki criterion [l3[ is decisive for entanglement of 
two-mode GS, and it can be given in terms of the quan- 
tity: 

S (V fi ) = hh + (1/4 - |/ 3 1) 2 ~h- 1/4(7! + 7 2 ), (3a) 

where 7i^ = detVi,2, h = detC, and I4 = 
tr IViZCZVzZC* Z] are invariants under local unitary 
operations, with Z = diag{l,— 1}. A GS is separable if 
and only if 

S(Vp)>0 7 (3b) 

the so-called Simon criterion. 

Entanglement in two-mode systems reminds Einstein- 
Podolsky-Rosen (EPR) discussion on completeness of 



quantum mechanics (l5| . Indeed, entanglement in this 
system is closely related to squeezing in "EPR-likeSS 
quadratures. Given any local operators Qi,Pi satisfying 
the commutation relations [Qi, Pi] = ih, for any arbitrary 
real number a ^ one can define the pair of EPR-like 
operators u — \a\Q\ + \Q2-, v = \ a \Pi — -Pi- Duan 
et al. (22l | have shown that if a state is separable, then 
(Am 2 ) + (Av 2 ) > a 2 + l/o 2 . Hence, if this pair of EPR 
operators is squeezed enough, i.e., if the sum of their 
variances violates the inequality, the state is entangled 
for sure. 

For GS, this criterion is necessary and sufficient to de- 
cide separability, which we shall call the Duan-Giedke- 
Cirac-Zoller (DGCZ) criterion. This is done representing 
the CM in a standard form, through local Gaussian oper- 
ations, so that the validity of the above inequality applied 
to the matrix in this form, with a determined by its co- 
efficients, implies that the state is separable. Restricting 
further to the set of symmetric states, it is sufficient to 
consider \a\ = 1, and the procedure amounts to finding, 
through local rotations and squeezing of the quadratures 
Xi,pi, which pair of EPR-like operators have the least 
value for the sum of their variances. 

We shall deal mostly with symmetric GS (i.e. I\ = I 2 ), 
and for these we shall use the Entanglement of Forma- 
tion as an entanglement quantifier, which has an explicit 
formula [23[ 

E F (p) = f(2^h + \h\- ^h + 2h\h\) (4) 

where f(x) = c+(x) log 2 c+(x) — c-(x) log 2 c_(x). with 
c±{x) = \{x- 1 l 2 ±x 1 ' 2 ) 2 . 



B. Two- mode entanglement from engineered 
reservoir 

To make the context clear, we now review the scheme 
for constructing the common squeezing reservoir between 
the modes . We also take the opportunity to introduce 
notation and to explicit the dependence of the final mas- 
ter equation on the several parameters of the setting. 

Consider two modes of electromagnetic field of a high 
quality microwave cavity, with frequencies u>\ and L02 ■ As 
before, a,; denotes the annihilation operator for mode i. 
The engineered reservoir is provided by a stream of atoms 
passing through the cavity. The atoms are first prepared 
in a specific superposition of two Rydberg states denoted 
\g) and |e), then they pass through the cavity, where they 
can interact with the two nondegenerate modes, while a 
classical field (injected externally in the setup of open 
cavities) saturates the dipole transition, pumping the two 
modes. The Hamiltonian that describes such setup is: 

H = huj Q a + a- + hn(e- iULt a + + e l " Lt a~) 

+ ^2 l [fk J j i a[ai + hg t (a,ia + + ajo-~)], (5) 
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where luq is the transition frequency between the atomic 
levels, gi are the coupling constants between the atom 
and the modes. a + and a~ are the atomic ladder opera- 
tors. The coupling with the external classical field, with 
strength O, is described by the time dependent part. For 
future reference, we set A = lul — ujq, the detuning be- 
tween the classical field and the atomic level. 

The authors explore different approximations and 
mode redefinitions. Here we will be concerned with the 
regime when (i) the atomic coupling with the classi- 
cal field is much stronger than with the cavity modes: 
1^1 ^ I Si I; (ii) defining d = \/A 2 + 40 2 , and choosing 
u>l obeying lul — cui = cu-2 — ujl = d, with the condition 
g = gi = g2- Under this regime the interaction Hamilto- 
nian can be approximated by: 

H int - -m b (bi#- + b\n + ), if A > 0, (6a) 
H mt ~ m b (bln- + b 2 n+), if A < 0, (6b) 

where n + and 7r~ are ladder operators for semiclassi- 
cal dressed states |+) = sm9\g) + cos#|e) and |— ) = 
cos 6 \g) — sin 9 \ e) , with tan 8 = 20/ (d — A) , Ob is related 
to the coupling constant between the atoms and the cav- 
ity modes, Ob = g^J (1 — /i) /(f + /i), where p, = tan 2 9 
\fi = (tan6>)~ 2 ] if |tan0| < 1 [|tan0| > 1], i.e., p, is 
determined by the classical field parameters. The new 
modes are 6 1(2 ) = S' t (r AI )a 1 ( 2) S'(r A1 ) = cosh |r AI |a 1(2 ) - 
sinh Ir^la^i) defined by the well-known two-mode 

squeezing operator: S(r^) = exp(r*aia,2 — r^a\a\), and 
r M is the squeezing parameter = arctanh^t. 

If after the interaction time the atoms are ignored, 
the Hamiltonian of Eq. ([5]) implies an open system ef- 
fective dynamics for the field modes. Eqs. ([6]) show that 
for A > (A < 0), the interaction reduces to a (anti-) 
Jaynes-Cummings between the classically dressed atom 
and mode &i(2)- F° r A > (A < 0), one simulates 
null temperature dissipation in mode &i(2) if one prepares 
atoms in the state |+) (|— )) and allow interaction for 
times, t, obeying ObT <C I. We shall call these type 1 (2) 
atoms. If a stream of type j atoms passes trough the cav- 
ity, one at a time, the field dynamics will be Markovian, 
given by a differential equation in Lindblad form (24) : 

■Jj; =Vj,Eng[p), (7a) 
where the effective engineered dissipator is given by |25| 
Vj,En g (p) - 2«,-(26,-pS] - btbjp ~ pb}bj), (7b) 

while Kj = (r a t ,j0 2 r 2 )/4 with r at j the atomic prepara- 
tion rate. 

We note that, at first sight, the large frequency sep- 
aration between the modes would forbid, in the regime 
considered here, the presence of combined terms between 
the modes, such as aia\ in the master equation, since 
they are fast oscillating compared to the total time scale 
of the experiment and even the interaction time of each 



atom [which is essential, in order to use the approxi- 
mated Hamiltonian ©]. Though, the crucial point is 

that the terms , the only combined ones present 

in the master equation, do not oscillate in the laser refer- 
ence frame. We will come back to this point later when 
discussing the state evolution. 



III. ENTANGLEMENT DYNAMICS UNDER 
ENGINEERED AND THERMAL RESERVOIRS 

If a random source defines together the type of the 
atom and the suitable DC electrical field in the cavity, 
the dissipator for the engineered reservoir will acquire the 
form: 

T>Eng(p) = 'Dl,En a {p) + ^2.Eng(p), (8) 

with Kj = (r QtJ 2 T 2 )/4, r at .j being the type j atoms 
flux. 

The entanglement behavior will get richer when con- 
sidered together with the natural dissipation and ther- 
mal noise on the modes, effects that, in usual experi- 
ments with microwave cavities, can be well described by 
a Lindblad equation with dissipator T>Nat given by: 

T^Natip) = M n ?« + l){2a,ipa\ - a\a,ip - pa\a,i) 

i 

+X l n Tt (2a\pa,i - a 4 a-/5 - paia]), 

(9) 

where denotes the number of thermal photons and 
Xi the decay rate for mode i. Since the number of ther- 
mal photons is approximately the same for both modes if 
\u)i — W2I we shall assume from now on = TVr. 

At last, the full dynamics of the system will be de- 
scribed by the Lindblad equation: 

-j t =V Nat (p)+V Eng (p), (10) 

with the dissipator given by Eqs. (J5J and ©. 

In an experimental setting one should actually con- 
sider the evolution in the laser reference frame. In this 
case, fast oscillations (with frequencies of the order d) 
between the modes will take place. As a consequence, 
in the course-grained scale necessary for the approxima- 
tions to be valid, some coherences of the density matrix 
will vanish, making the theoretical analysis more compli- 
cated. Nevertheless, the above equation will still be valid 
for some initial states, i.e., they are not affected by these 
oscillations. Hence, in the rest of this section we shall ex- 
plore the theoretical properties of Eq. (fTTjj) , considering 
several initial states, but later, in the experimental pro- 
posal, we must guarantee the use of such robust states 
against this coarse-graining effect. 
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A. Symmetric engineered reservoir 

Now we pass to the study of entanglement dynamics 
under Eq. (fTU)l . Given the specific form of the dissipa- 
tors, Gaussianity is preserved, so we are safe to restrict 
ourselves to GS. Let us suppose first that both types of 
atoms enter in the cavity with equal probability, so that 
Ki = K 2 = and also, for the sake of simplicity, that 
Ai = A2 = A. The equations of motion for the second 
momenta in this case become: 



n 
m 
m 
m 



-2{k + X)rij + 2k\B\ 2 ■ 
-2(k + X)mj, 
-2(« + X)m c + 2kAB* 
-2(k + X)m s , 



2Xn T , 



(11a) 
(lib) 
(11c) 
(lid) 



with A = cosh(r), B = e 1 ^ sinh(r), where r is the squeez- 
ing parameter and (j> the squeezing angle in the definition 
of modes bj (r M = re 1 ^). This gives a relaxing dynamics 
with the asymptotic CM: 



Vi 



( n U + 5 
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m cJ 


n hf + \ 
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n 2.f + 5 











n 2 j + 



(12) 



\B\ 2 +n T R 



AB* 
l+R> 



and the 



where mj = n 2 j = — j+t^- 
ratio R = X/k was introduced. 

Simon criterion can now be applied to determine 
whether such CM represents entangled^ or separable 



states. In the same reasoning as in Ref. |H, S (V nt ) > 



implies a deep separable asymptotic state, i.e., a state 
belonging to the interior of the separable states set. This 
can be seen noting that S is a continuous function, so 
there must exist a "ball" of GS around the asymptotic one 
such that S is strictly positive, that is, a ball of separable 
GS. Hence, this translates dynamically as sudden death 
of entanglement (SDE), because we can say for sure that, 
for every initial GS, there will be a time T such that the 
entanglement will be zero for t > T. In principle, though, 
it could undergo cycles before it vanishes for good, or an 
initially separable state can acquire some entanglement 
and (necessarily) lose it after some time. These "non- 
asymptotic" behaviors will be discussed in more detail 
on Sec. llnBl 

On the other hand, states satisfying S (V Pf ) < repre- 
sent a situation of (asymptotic) persistent entanglement 
(PE) . Note that entanglement can be created by the com- 
mon reservoir, since this is exactly the idea of the Ref. (gj 
proposal. But again, the intermediate dynamics can ex- 
hibit richer features: for instance, an initially entangled 
state can lose all entanglement at finite time but (neces- 
sarily) recover it at later time, as will be exemplified at 
Sec. IITTB1 

The exceptional situation is given by S (V Pf ) = 0, 
when each initially entangled state can show one of two 
fates: sudden death of entanglement or asymptotic death 
of entanglement, depending on the initial state. Contrary 



to the former two situations, this one requires the knowl- 
edge of the whole dynamics in order to determine which 
fate will occur. 

The situation here studied allows only one asymptotic 
state for each set of fixed parameters, that is why one can 
not see infinite cycles of birth and death as in Ref. |f| , 
which would appear for dynamics with asymptotic peri- 
odic orbits, instead of asymptotic states. 

In Fig. [TJ we plot the regions in the parameter space 
Rxtit where the asymptotic state is separable or entan- 
gled, defining the fate of entanglement. The boundary 
curve separating the two regions in the diagram has the 
simple form: 



„2r 



1 



tit 



2R 



(13) 



The physical interpretation is simple and meaningful. 
Given the engineered reservoir, for any given positive 
coupling ratio, R, there is a positive temperature, tit, 
obeying (fT3|) such that below this temperature, the 
asymptotic state is entangled, due to the common reser- 
voir, while for temperatures above that critical value, the 
asymptotic state is separable, when (local) thermal noise 
prevails. 

For any set of reservoir parameters, we applied the 
method of Duan, Giedke, Cirac, and Zoller [22j and 
found the EPR-like operators with least sum of vari- 
ances, which gives Xi jt j, — ^2,0 and + P2,4>i where 
(X i ^,Pi t< j > ) T = 1l2<f,(xi,Pi) T , 7^20 is the matrix repre- 
senting a rotation in the plane trough an angle 2(f>, <f> be- 
ing the squeezing angle defined by the reservoir. This is 
expected, since the engineered reservoir tries to lead the 
initial state to a usual two-mode squeezed state, known 
to be squeezed in these quadratures. The natural reser- 
voir, on the other hand, enlarges their spreading. So 
the final decision of whether the asymptotic state will be 
separable or entangled will depend on this competition 
between the reservoirs: one trying to squeeze the collec- 
tive quadratures, the other trying to spread them. 

Since the asymptotic states here obtained are also sym- 
metric, we can apply Eq. (@| to calculate their entangle- 
ment of formation. In Fig. ® we show their values as 
functions of R for some fixed values of nx- For null tem- 
perature, when there is no thermal photon, the entangle- 
ment is positive for any rate. However, for each positive 
temperature, there is a maximal rate above which the 
state is deep separable (as the coupling to the natural 
reservoir grows, thermal effects are more sensible) corre- 
sponding to the regime where the system exhibit SDE. 



B. Non-asymptotic dynamics 

We will discuss here in more detail the system inter- 
mediate dynamics. This task is facilitated recognizing, 
by Eqs. (fTTj) . that the dynamics will always describe a 
straight line in the space of the parameters of the CM 
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FIG. 1. Representation of the asymptotic states determined 
by reservoir parameters nr and R with respect to their en- 
tanglement, for r — 1 and arbitrary (j>. Below the curve, they 
are entangled, and above, separable. This also determines 
the asymptotic properties for the entanglement dynamics of 
any initial state: for points above the curve, every initially 
entangled state exhibits SDE; for points below the curve, ev- 
ery initial state has some entanglement for large times; and 
for points exactly in the curve, entanglement can die at finite 
time or asymptotically. 



senting the values of n, m c such that the corresponding 
CM is nonphysical, separable or entangled. The frontier 
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FIG. 3. Representation of the CM's with m = n% = n, m c > 
and all other elements equal to zero. The region where these 
parameters represents physical CM's is divided by a straight 
line, above which all states are entangled, while below they 
are separable. Both regions are convex in the parameters. 
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FIG. 2. Entanglement of the asymptotic state as function of 
the ratio R — A/ft between the reservoir's coupling constants 
for some values of the number of thermal photons, tit, with 
r — 1 and arbitrary <f>. While the entanglement is always 
positive for tit = 0, if tit > there is always some value of 
R above which the entanglement is zero. 



(ni,ri2, etc.)[? ], for every initial GS, and that the set 
of separable states is convex for these parameters also. 
In other words, the trajectories of the CM are given 
by straight line segments exponentially approaching the 
asymptotic CM, Vp. . From now on we set </> = 0, for 
simplicity in the analysis. 

We begin with the situation where the reservoir param- 
eters satisfy Eq. (|13p so that the corresponding asymp- 
totic state rests in the border between the separable and 
entangled sets. If we restrict our attention to initial states 
with CM, where n\ = ri2 = n, m c is real and all other 
elements are null, Eqs. ([lip keep the dynamics inside 
this same subset. In Fig. [3J we plot a diagram repre- 



between separable and entangled states is given by the 
straight line n = m c , so the entangled region is convex 
in these parameters, and by our choice an asymptotic 
state is denned by the dynamics somewhere in this line. 
A simple picture of the dynamics can be given for these 
states, using the fact that the trajectory in the n x m c 
space will be a straight line, from the values (no,m C fi) 
of the initial state to the asymptotic values (n/,m Cj /), 
defined by reservoir parameters r and R [with n-r given 
by Eq. (JT3J)] . If these parameters are such that (no, m c ,o) 
represents an entangled state, since (rif,m c _f) belongs 
to the line n — m c , the entire trajectory will be on the 
entangled-state region, by convexity, i.e. entanglement 
will vanish asymptotically. To illustrate this, we plot 
in Fig. U] (dashed line) the evolution of the S function, 
parametrized by p(t) = 1 — cxp [— 2(n + X)t], for an initial 
state with uq = 1 and m Cj o = 1.0125 and an asymptotic 
state with nf ~ m c j = 1. Its value is initially negative, 
since the initial state is entangled, and remains negative 
for all times. 

For examples of entanglement sudden death, consider 
the set of states with n\ = n-2 = n, m c = 1 and 
mi = i~n-2 = ni a real number. In Fig. [5l we exhibit a 
diagram analogous to the one in Fig. [31 and we see that 
the subset of entangled states is not convex anymore. If 
the asymptotic state also has m c = 1, the trajectory of 
the state CM can be described in this diagram and will 
be again a straight line. Since the entangled region is not 
convex, we may have initially entangled states that will 
lose all entanglement at finite time, even if the asymp- 
totic state is in the frontier between the regions. For 
instance, we can take parameters for the reservoir such 
that n/ = m c j = 1 and an initial CM that has elements 
no = 1.2, toq = 0.5, and m c .o = 1- In Fig. U (continuous 
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FIG. 4. Evolution of the Simon function S, parametrized by 
p(t) for three distinct situations. The continuous line cor- 
responds to an initially entangled state, which enters in the 
separable region at finite time and converges to an asymp- 
totic state in the frontier. The dotted line is for the same 
initial state, but the reservoir is such that its corresponding 
asymptotic state is entangled, with the system losing all its 
entanglement to recover it later. Finally, the dashed line cor- 
responds to a situation of asymptotic death of entanglement. 
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FIG. 6. Evolution, parametrized by p, for the sum of vari- 
ances of pairs of EPR quadratures, subtracted by 2, so that a 
negative value indicates an entangled state. The continuous, 
dotted, and dashed curves are, respectively, for the optimum 
pairs of the final state, the initial state, and the state in each 
instant of time. We see that the common reservoir enlarges 
the optimum pair of the initial state to squeeze the pair of 
the asymptotic state (which is initially large). And there is 
a period of time where the quadratures with the least sum of 
variances are not squeezed enough to entangle the states. 
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FIG. 5. Representation of the CM's with ?ii = n% = n, mi — 
m,2 = m > 0, m c — 1 and all other elements equal to zero. 
The region where these parameters represents physical CM's 
is divided by the dashed curve, above which all states are 
entangled, while below they are separable. Now the entangled 
region is not convex in the parameters. 

line), we plot, as before, the function S of the evolved 
state. We see that it is initially negative, representing an 
entangled state, it became positive at finite time, so the 
state has no entanglement, and it remains strictly posi- 
tive until it vanishes for p = 1 (or t — > oo), because the 
state converges to a point in the frontier. 

If we perturb a little the reservoir parameters so that 
this asymptotic state becomes entangled, and taking the 
same initial state, the curve in the figure will be slightly 
distorted but now will cross the p axis twice, meaning 
that the state loses all entanglement suddenly but re- 
covers it a later time, remaining entangled for the rest 
of the dynamics, a possibility mentioned before (see the 
dotted line of Fig. |H with the parameters n/ = 0.95 and 



m c j = 1 for the asymptotic state). This is a consequence 
of the fact that the set of entangled states is not convex 
on these parameters. Cycles of birth and death are not 
allowed, since a straight line can cross the convex set of 
separable states only once. 

This feature can be understood, from a more physi- 
cal point of view, considering the entanglement from the 
perspective of squeezed EPR-like quadratures. Again, 
the EPR-like operators with least sum of variances val- 
ues (i.e., the optimum pair for this state) can be found 
and read Xiy — X^y and P\y + Pa,' with Xiy = r'&i, 
Pt,r> = Pi/r' and r' = y/{n - m - 1/2) /(n + m - 1/2). 
But these are not the quadratures "chosen by" the com- 
mon reservoir to squeeze (being those with r' = 1, since 
we had set cj> = 0). So, on its way to squeezing its fa- 
vorites, it enlarges the ones from the initial state, such 
that there is a period of time when no pair of EPR-like 
quadratures at all are squeezed enough to entangle the 
modes (see Fig. \Q. 

It remains then to explore the intermediate dynamics 
for the situation where the asymptotic state is in the in- 
terior of the separables, but this is also easily inferred 
from the straight line trajectories exhibited by the sys- 
tem: every initial separable state will remain separable 
for all times, and every initial entangled state will lose 
its entanglement in finite time. 



C. Asymmetric engineered reservoir 

We have also considered the other extreme case where 
only one type of atom enters in the cavity, say of type 
1, so that ki > and Hi — (but again Ai = A2 = A). 
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Here, the second momenta equations of motion read: 

rii = -2(A 2 n + X)m + (ABn)m c + 

+(AB*K)m* c + 2Xn T , (14a) 
ri 2 = -2(X - |B| 2 K)n! - tti c {ABk) ~ 

-(AB*n)m* c + 2Xn T + 2n\B\ 2 , (14b) 

mi = -2(A 2 k + X)mi - (2AB*K)m s , (14c) 

m 2 = -2(X - \B\ 2 n)m 2 + (2AB*K)m*, (14d) 
m c = — (k + 2X)m c — (AB* n)ri\ + 

+ (AB*n)n 2 + AB* k, (14e) 
?ii s = -(re + 2X)m s + (ABre)mi - (AB* n)m 2 . (14f) 

Despite the asymmetry of the master equation here, 
the asymptotic states have the same qualitative behavior 
as in the previous case: for null temperature it will be 
entangled for any value of R = rei/A but for finite tem- 
perature it is separable for high enough R. In Fig. [7] we 
plot the diagram analogous to the one in Fig.[T] determin- 
ing the dynamical phases for the entanglement dynamics, 
also separated by a curve with a similar shape. 



2.0 r 



1.0 - 




0.0 h 

1 2 3 4 5 

R 

FIG. 7. Phase diagram for the asymmetric reservoir, exhibit- 
ing the same qualitative behavior as in the symmetric case. 

The asymptotic states entanglement, using now the 
log-negativity for asymmetric GS pjj], are also shown 
(Fig. [8} for distinct values of tit, as function of R, and 
the same interpretation applies here. 

Analogous results are valid also when k 2 > but 
Ki =/= k 2 . So, as far as asymptotic issues are concerned, 
the system presents the same general behavior for both 
symmetric and asymmetric reservoirs. However, the non- 
asymptotic analysis made before for the symmetric reser- 
voir becomes much more involved for the asymmetric one, 
and will be subject of future work. 



IV. EXPERIMENTAL PROPOSAL 

The dynamics studied here is simpler than that in 
Ref. |f| (e.g. being Markovian), yet it is much more con- 
trollable and with clear experimental motivation. In this 




0.0 0.5 1.0 1.5 2.0 

R 



FIG. 8. Log-negativity of the asymptotic state for the values 
r = 1, and some values of ixt- 

section, we discuss in more detail how to test these pre- 
dictions in the laboratory. 

Considering the state evolution in the laser reference 
frame, Eqs. (fTT|) will be replaced by 

rij = -2(« + X)nj + 2n\B{ 2 + 2Xn T , (15a) 

rhj = -2(k + X)mj + (-l) j+1 2idm j , (15b) 

m c = -2(k+ X)m c + 2kAB*, (15c) 

m s = — 2(k + X)m s + 2idm s , (15d) 

where only the equations for m,j and m s have changed. 
If we consider the initial GS where rrij = m s = 0, as 
well as null amplitudes, the evolution will be the same as 
that for the interaction picture, i.e., given by Eq. (|10[) . 
and the coarse graining in time will not affect these 
states throughout the evolution. For other initial states, 
the coarse graining will be a source of decoherence, and 
the state may even lose the Gaussian character along 
the evolution (like a one-mode coherent state, with non- 
vanishing amplitude, averaged over its phase). 

Our proposal starts then from the preparation of a 
highly entangled state using the common symmetric 
reservoir Q for a total time much shorter than A" 1 , and 
setting R <C 1 in order to make thermal effects negligi- 
ble. After the preparation time, the parameter R is set 
to values of the order R ^ 1, by changing the flux of 
atoms through the cavity [? ]. The system is allowed 
to evolve until some time t (of the order A -1 ) and the 
entanglement evolution can be studied, for example, by 
using the method proposed in Ref. (26| , to reconstruct the 
two- mode state. Naturally, one could also try to measure 
some entanglement witness to simplify the experimental 
procedure, by using the previous knowledge of the field to 
reduce the amount of experimental data, but still char- 
acterize its entanglement. 

For this setting, i.e., assuming an initial two-mode 
squeezed state, the times where entanglement sudden 
death (ESD) takes place can be obtained explicitly and 
reads: 

R 

Xt E SD= „ n , D , ln(l -Pesd)^ 1 , (16) 
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where 

(1 + R)(B 2 -AB) 
PESD (1 + R)(B 2 — AB) — B 2 — utR + AB ' 1 ' 

and the parameter R refers to the one in the second 
stage of the procedure, after the initial state preparation. 
These times, in units of A -1 , are exhibited as function of 
R in Fig. [9j for some values of temperature. Naturally, 
they are infinite for small R, in the region where per- 
sistent entanglement takes place. After the threshold, it 
falls abruptly and stabilizes in values close to unity, so 
the entanglement will typically survive long enough for 
its (not so sudden) death to be observed. It may also 
be interesting to study thermal effects like the behavior 
of ^f- before the sudden death, which is an indicator of 
the incidence of the dynamical trajectory of p at the set 
of separable states, by controlling the temperature. 



5r 




2 4 6 8 10 12 14 

R 

FIG. 9. Sudden death time as function of R, for some values of 
fix • It is infinite for small R, corresponding to the persistent 
entanglement region, but after the threshold, it rapidly falls 
to a value near to unity. 



V. ROBUSTNESS OF THE SCHEME FOR 
GENERATING ENTANGLED GS IN A CAVITY 

Since entanglement in GS is related to squeezing in 
the proper EPR quadratures, while the natural ther- 
mal effect is to spread Wigner functions, it is in order 
to ask about the sensibility of the preparation scheme 
with natural dissipation, considering the regime fl< 1 
and tit <C 1. Actually we have analyzed a slightly dif- 
ferent procedure than the original one by Pielawa et al. 
[f| . There they propose to empty the modes bj at turns, 
passing a stream of atoms of type 1 , say, until mode b\ is 
sufficiently "washed", then repeating the procedure with 
type 2 atoms. Here we consider that both types of atoms 
can pass through the cavity, not at the same time, but 
with equal probability, so that we can apply the equa- 
tions for the symmetric reservoir to compute the system 
evolution. 

Supposing that the two original modes, hi, are initially 
empty (and remembering that the number of photons 



can be lower than the number of thermal photons if one 
passes a suitable stream of atoms to "steal" photons from 
the mode) , the evolution of the covariance matrix of the 
system will be then given with the following non-zero 

entries: n\ = ni = jj^f-^ pft) and m c = j^p(t), 
where again p(t) = (1 — exp [— 2(k + A)i]). Assuming 
also that the duration of experiment t is large enough so 
that p»l (which can be achieved with f>K _1 but still 
t <C A -1 ) or, in other words, that the system is essen- 
tially in the asymptotic state of the procedure, we plot 
in Fig. [TU] the entanglement of formation as a function 
of R, for tit = 0.05 (value attained in the recent ex- 
periment reported in Ref. (2Tj ) for distinct values of r. 
The graphs are normalized by the value of the entangle- 



e f (R)_ 

E F (R=0) 




0.2 - 



0.0 0.1 0.2 0.3 0.4 0.5 

R 

FIG. 10. Entanglement of the asymptotic states as function of 
the ratio R — X/k between the reservoir decay rates, for some 
values of squeezing parameter r, with nx = 0.05, normalized 
by its value at R = 0. The continuous line is for r = 1, dotted 
for r — 1.5, dashed for r — 2 and dotted-dashed for r = 2.5. 
The entanglement is somewhat sensitive to dissipation, and 
is more sensitive the greater the squeezing parameter. 

ment of formation with R — 0, which would be obtained 
if there were no dissipation. We see that the entangle- 
ment is somewhat sensitive to the dissipation, even at 
such a low temperature, and can be lowered by half, for 
R = 0.1, i.e. , when the engineered reservoir rate is ten 
times greater then the dissipation rate. Also, the greater 
the squeezing parameter, the sensitive entanglement be- 
comes to dissipation, as expected. Though, if k is above 
two orders of magnitude greater than A, the entangle- 
ment does not appear to be significantly changed. For a 
squeezing parameter of r = 1, a value of R ~ 10~ 2 and 
a waiting time of about 3k _1 would suffice to obtain an 
amount of entanglement greater than 90% of the pure 
ideal state (with R = and infinite waiting time). 



VI. CONCLUSION 

An engineered reservoir can be used to create entan- 
gled states of two field modes in a cavity jgj. We here 
show that it can also be used, together with the natu- 
ral thermal reservoir, to depict asymptotic entanglement 
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phase diagrams Taking advantage of the fact that 
the specific reservoirs used preserve Gaussianity, we could 
make all the discussion of entanglement in terms of co- 
variance matrices. We showed that only two behaviors 
are allowed: asymptotic entanglement or entanglement 
"sudden death". The line between both phases allows the 
coexistence of both behaviors. 

Besides asymptotic studies, the symmetric engineered 
reservoir was also studied in detail, with examples of all 
possible entanglement fates exhibited. 

This study can be viewed as an experimental proposal 
for drawing entanglement phase diagrams and for moni- 
toring the "transition". Another byproduct is the study 



of the robustness of the proposal for entanglement gen- 
eration with respect to the natural thermal environment. 
An interesting question raised is what would be the best 
setup for such experimental drawing of asymptotic en- 
tanglement phase diagrams. 
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